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MOLECULAR INSIGHTS INTO THE BIOGEOGRAPHY AND 
SPECIES STATUS OF NEW ZEALAND’S ENDEMIC 
LATRODECTUS SPIDER SPECIES; L. KATIPO AND 
E ATRITUS (ARANEAE, THERIDIIDAE) 
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ABSTRACT. New Zealand’s endemic sand dune Latrodectus widow spider species, L. katipo and L. 
atritus , possess behavioral and physiological attributes likely to promote dispersal over large distances. 
Morphological, physiological and behavioral similarities between L. katipo and L. hasselti , an Australian 
endemic, suggest gene flow may occur across the Tasman Sea. In this study we examine intraspecific and 
interspecific genetic relationships within the ND1 gene region between L. katipo , L. atritus , L . hasselti 
and L. hesperus to assess whether the genetic evidence supports current taxonomic species designations. 
We found low interspecific pairwise distances among L. katipo and L. atritus populations, suggesting 
either introgression, incomplete lineage sorting, or that the current taxonomic distinction between the two 
species may be invalid. Parsimony and maximum likelihood analyses were inconclusive as to the rela¬ 
tionships between the New Zealand Latrodectus species and the Australian L. hasselti. Low pairwise 
distances between L. hasselti and the New Zealand widow fauna indicated that L. katipo and L. atritus 
were not present in New Zealand before the fragmentation of Gondwana. 

Keywords: Latrodectus , New Zealand, Australia, dispersal, molecular phylogenetics 


New Zealand’s endemic Latrodectus Wal- 
ckenaer 1805 fauna is considered to comprise 
two endemic species, L. katipo Powell 1870 
and L. atritus Urquhart 1890 (Forster & For¬ 
ster 1999). Latrodectus atritus was originally 
described as a subspecies of L. katipo (Ur¬ 
quhart 1890) and was proposed as a subspe¬ 
cies of L. hasselti Thorell 1870 by Parrott 
(1948), and a junior synonym of L. mactans 
(Fabricius 1775) by Levi (1959). McCutcheon 
(1976) and Forster & Kingsford (1983) argued 
that L. atritus is a separate species to L. katipo 
and Forster (1995) elevated L. atritus to spe¬ 
cies but did not provide any taxonomic justi¬ 
fication. The only reported morphological dif¬ 
ference between L. atritus and L. katipo is 
coloration, which is usually unreliable for sep¬ 
arating spider species but can be useful in sep¬ 
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arating Latrodectus species (McCrone & Levi 
1964). Latrodectus atritus females do not 
have the red median stripe on the dorsal sur¬ 
face of the abdomen that L. katipo has (Mc¬ 
Cutcheon 1976; Forster & Kingsford 1983; 
Forster & Forster 1999) and the males of the 
two species have slight differences in color 
(Forster & Kingsford 1983). Forster & Kings¬ 
ford (1983) also reported differences between 
the species in the time it took for spiderlings 
to emerge from the eggsac. Forster & Forster 
(1999) noted that L. atritus eggs and spider- 
lings need higher temperatures than L. katipo 
and they also stated that “laboratory studies 
show that they do not generally crossmate but 
when they do, the eggs are infertile” but no 
data was included in the publication to support 
this. Both New Zealand species inhabit coastal 
dune systems and commonly build webs in 
low growing plants and driftwood or flotsam. 
Although the niches occupied by L. katipo and 
L. atritus are similar, their known distributions 
are distinct. Latrodectus katipo inhabits coast¬ 
al dunes in the northern half of South Island 
and the southern half of North Island, whereas 
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L. atritus inhabits coastal dunes in the north¬ 
ern half of North Island. There is an overlap 
in the species 5 distributions in Taranaki on the 
west coast of North Island (Forster & Forster 
1999) and in Hawk.es 5 Bay on the east coast 
North Island. 

Discriminating between Latrodectus spe¬ 
cies using morphology has always been prob¬ 
lematic (Levi 1983) and the distinction be¬ 
tween L, katipo and the Australian L, hasselti 
is minimal. There is no difference between the 
male pedipalp and female external and inter¬ 
nal genitalia of the two species (Parrott 1948; 
Levi 1959) but that is often the case in Lat¬ 
rodectus (Levi 1983). Forster & Forster 
(1999) state that the most definitive morpho¬ 
logical character that separates the two species 
is the dense covering of short, fine hairs on 
the body of L katipo compared to the long 
fine hairs and stouter short hairs on L. hasselti . 
Latrodectus hasselti females are usually larg¬ 
er, their webs are stronger and usually yellow¬ 
ish in color (Forster & Forster 1999) and this 
species performs a stereotyped behavior of 
sexual cannibalism (Forster 1992). Laboratory 
studies examining interspecific interactions 
have shown that L. hasselti females will not 
mate with L. katipo males (Forster 1992, 
1995). However, in their phylogenetic study 
of Latrodectus , Garb et al. (2004) found that 
L katipo and L. hasselti were closely related 
(4.9% uncorrected divergence in a 428 bp sec¬ 
tion of the mitochondrial gene cytochrome ox¬ 
idase subunit I (CGI)). 

Recent research demonstrates that L. katipo 
and L. atritus are probably good dispersers. 
Spiderlings of both species are able to dis¬ 
perse by ballooning and adult females may be 
able to disperse on driftwood at sea under op¬ 
timal conditions (Griffiths 2002). It is uncer¬ 
tain how 7 far L. katipo , L, atritus , or L. hasselti 
could disperse by these means, but other bal¬ 
looning spider species have been recorded 
landing on ships up to 300 km from land 
(Gertsch 1979), indicating that ballooning spi¬ 
ders may travel substantial distances. Further¬ 
more, spiders recorded from driftwood and 
flotsam at sea, suggest that w T ater-bome spi¬ 
ders could also disperse over large distances 
(Foelix 1996) and adult L. hasselti can survive 
more than 300 days without food (Forster & 
Kavale 1989). This evidence may explain why 
L. katipo and L, atritus distributions span nu¬ 
merous geographic barriers such as headlands, 


estuaries, rivers and areas of open sea (<30 
km) and may account for morphological, 
physiological and molecular similarities be¬ 
tween the New Zealand Latrodectus fauna and 
L. hasselti , which is considered endemic to 
Australia (Parrott 1948; Levi 1959; Forster & 
Kingsford 1983; Garb et aL 2004). As yet, 
however, the influence dispersal may have had 
on the biogeography of L. katipo and L. atri¬ 
tus has not been investigated. Moreover, be¬ 
havioral, morphological, physiological and 
molecular similarities between L. katipo and 
the Australian endemic, L. hasselti , have not 
been adequately explained. 

In this paper, we examine intraspecific and 
interspecific genetic relationships within the 
NADH dehydrogenase subunit 1 (ND1) mi¬ 
tochondrial gene region between L. katipo , L, 
atritus , L. hasselti and L. hesperus Chamber¬ 
lin & I vie 1935 (all part of the strongly sup¬ 
ported “mactans clade” In Garb et al. 2004) 
to assess the degree of separation between the 
New Zealand and Australian Latrodectus spe¬ 
cies and whether genetic evidence supports 
current taxonomic species designations. The 
ND1 gene region was chosen because it is fast 
evolving and has been successfully used to ex¬ 
amine genetic differences between spider spe¬ 
cies and populations (Hedin 1997a, 1997b; 
Masta 2000; Johannesen et al. 2002; Maddi- 
son & Hedin 2003; Masta & Maddison 2002; 
Vink & Paterson 2003). 

METHODS 

Adult female L. katipo and L. atritus were 
collected from eight sites around New Zea¬ 
land (Fig. 1) and w ? ere stored in 95-100% 
EtQH at “”80 °C to maintain high quality 
DMA. Voucher specimens are stored at the 
Ecology and Entomology Group, Research 
Collection, Lincoln University, New Zealand. 
Specimens were collected from sites that were 
selected throughout the distributions of both 
species. One specimen per population of L, 
katipo was collected from Kaitorete Spit 
(43°50'S, 172°31'E) and Waikuku Beach 
(43°17'S, 172°43'E), Canterbury, from Fare¬ 
well Spit (40°30'S, 172°48'E), Golden Bay 
and from Flat Point (41°28 ; S, 175°37'E) and 
Herbertville (4G°29'S, 176°37'E) on the east 
coast of the low 7 er North Island (Fig. 1). One 
specimen per population of L, atritus was col¬ 
lected from Houpoto (37°58'S, 177°33'E), 
Rarawa (34°44'S, 173°05'E) and Opoutere 
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(37°24'S, 175°56'E) in the upper North Island 
(Fig. 1). Latrodectus hasselti were collected 
from Myalup, Western Australia (33°06 ; S, 
115°41'E) and Brisbane, Queensland 
(27°27'5, 153°02). A specimen of L. hesperus , 
intercepted by the New Zealand Ministry of 
Agriculture and Forestry on table grapes from 
California, USA was used as an outgroup. Al¬ 
though L . hesperus is common throughout 
western North America (Chamberlin & Ivie 
1935; Levi 1983), the identification of the L. 
hesperus specimen is tentative as a second, 
undescribed Latrodectus species is reported to 
be present in California (see Levi 1983) and 
the taxonomic differences between the two 
species are unknown. The entire front leg and 
hind leg from one side of each specimen were 
removed and washed in sterile deionized, dis¬ 
tilled water to remove excess alcohol. Geno¬ 
mic DNA was extracted from samples using 
a proteinase-K digestion and high salt precip¬ 
itation method (White et aL 1990). The DNA 
was suspended in 1:20 TE (1 OrriM Tris, 1 rriM 
EDTA, pH 8.0). 

The first half (—420 bp) of the mitochon¬ 
drial NB1 gene region was amplified from di¬ 
luted genomics in 25 juil PCR reactions using 
the primers Nl-J-12261 and TL1-N-12718 
(Hedin 1997a). Each 25 fxl reaction contained 
1X Taq buffer, 1 mM dNTPs, 2 pM MgCl 2 , 
0.4jxM of each primer, l .25 units Taq DNA 
polymerase (Roche), and 1 pi diluted genomic 
DNA. Amplification took place in a Gene- 
Amp® 2400 Thermocycler and included an 
initial denaturation of 4 min. at 94 °C fol¬ 
lowed by 40 cycles of 40 s at 94 °C, 40 s at 
45 °C, 40 s at 72 °C and a final extension of 
5 min. at 72 °C. The resulting PCR product 
was purified by precipitation with 50 pJ of 
isopropanol and 25 pi NH 4 Ac (2.5M) to re¬ 
move excess salts and primers. Purified ds- 
DNA samples were washed in 70% EtOH and 
suspended in 6 pi of sterile deionized, dis¬ 
tilled water. All dsDNA samples were subse¬ 
quently sent to the Waikato DNA Sequencing 
Facility where they were sequenced in both 
directions. 

DNA sequences were aligned against a 
complementary-strand sequence in DNA- 
MAN (version 4.02), and checked against 
hard copy chromatograms by eye. Corrections 
were made where necessary. The possibility 
of pseudogenes and polymerase errors were 
eliminated by the translation of the sequences 


to amino acids and no stop codons or frame- 
shifts were found, A multiple alignment of all 
sequences was compiled in CLUSTALX 
(Thompson et al. 1997) and imported into 
PAUP* 4 0b 10 (Swofford 2002) for analysis. 

Data were analyzed as unordered characters 
using parsimony with the exhaustive option 
selected. Bootstrap values (Felsenstein 1985) 
for monophyletic groups were calculated us¬ 
ing the branch and bound search option in 
PAUP*. Model test version 3.06 (Posada & 
Crandall 1998) was used to select the maxi¬ 
mum likelihood parameters and the HKY + 
F model (Hasegawa et al. 1985) was used to 
estimate the maximum likelihood tree. The 
branch and bound option was selected in 
PAUP* for the maximum likelihood analysis 
and branches were collapsed creating polyto- 
mies if the branch length was < le—08. Boot¬ 
strap values for the maximum likelihood tree 
were calculated using a heuristic search 
(10000 replicates). Base frequency calcula¬ 
tions, transition/transversion ratios, number of 
variable sites and the conversion of nucleo¬ 
tides to amino acids were conducted using 
MEGA version 2A (Kumar et al 2001). 

In addition to the molecular work, 22 spec¬ 
imens of L, katipo and L. atritus from collec¬ 
tions at the Museum of New Zealand, Otago 
Museum, Auckland Museum and Lincoln 
University Entomology Research Museum 
were examined for differences in male and fe¬ 
male genitalia. 

RESULTS 

The nucleotide composition was G (gua¬ 
nine) depauperate (29% A, 22% C, 10% G, 
39% T), which is similar to that of the spider 
family Nesticidae (Hedin 1997a), a sister fam¬ 
ily of Theridiidae (Griswold et al. 1998). Se¬ 
quence data were deposited in GenBank (Ben¬ 
son et aL 2002), accession numbers 
AY383604-AY383614. 

The largest interspecific pairwise distance 
between L. katipo , L. atritus and L. hasselti 
was 1.95 %, whereas the smallest pairwise 
distance between the Australasian specimens 
and L. hesperus was 21.41 % (Table 2). In 
contrast, the largest intraspecific pairwise dis¬ 
tance between L. katipo specimens was 0.97 
%, which was the same as the largest pairwise 
distance between L. katipo and L. atritus 0.97 
% and greater than the largest pairwise dis- 
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Figure 1.-—Sites where L. katipo (O —unfilled circle) and L. atritus (#—filled circle) samples were 
collected for molecular analysis. 


tance between L, atritus specimens 0.73 %, or 
L. hasselti specimens 0.00 % (Table 2). 

Three haplotypes of L. katipo occurred 
among the specimens sampled. The L. katipo 


specimens from Hcrbertville and Flat Point 
had identical ND1 sequences as did the spec¬ 
imens from Farewell Spit and Wiakuku. There 
were three haplotypes of L. atritus and all 
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Table 1.—Base frequencies, transition/transversion ratios and number of variable sites at each codon 
position. 



Position 1 

Position 2 

Position 3 

Base frequencies % 

T 32.8, C 18.4, 

T 44.9, C 21.3, 

T 40.1, C 26.4, 


A 36.0, G 12.7 

A 21.9, G 11.9 

A 28.7, G 4.7 

Transition/transversion ratio % 

3.5 

1.3 

1.6 

Variable sites including L. hesperus 

16 

7 

68 

Variable sites excluding L. hesperus 

1 

0 

10 


three specimens had different sequences. The 
ND1 sequences of both J L hasselti specimens 
were identical. Ninety-one nucleotides varied 
among the four Latrodectus species but only 
11 sites varied between the L katipo , L. atri¬ 
tus and L. hasselti specimens. Table 1 lists the 
base frequencies by base position, transition/ 
transversion ratios and number of variable 
sites at each position. 

Fourteen of the 137 amino acids coded for 
by the ND1 sequence were variable. Only two 
were variable within L, katipo , L. atritus and 
L. hasselti . There were three Australasian hap- 
lotypes. The two L. hasselti specimens formed 
one haplotype. Latrodectus katipo from Flat 
Point, Herbertville and Kiatorete Spit had 
identical amino acids. The third amino acid 
haplotype consisted of L. katipo from Fare¬ 
well Spit and Waikuku and L. atritus from 
Houpoto, Rarawa and Opoutere. 

Parsimony analysis yielded 12 equally par¬ 
simonious trees, 97 steps long with a consis¬ 
tency index, excluding uninformative charac¬ 
ters, of 0.714 and a retention index of 0.714. 
Of the 411 characters included, 91 were var¬ 
iable of which four were parsimony informa¬ 
tive. There was no consensus among the 12 
trees. Maximum likelihood analysis yielded 1 
tree with a score of —in 883.451. The likeli¬ 
hood tree (Fig. 2) had an identical topology 
to one of the 12 most parsimonious trees. 
There was no bootstrap support for any of the 
clades in the parsimony analysis and only 
very weak bootstrap support in the maximum 
likelihood trees (not shown on Fig. 2); 59% 
for a clade containing Queensland and West¬ 
ern Australia, Farewell Spit and Waikuku, 
Houpoto, Gpoutere, Rarawa, and Kiatorete 
Spit and 52% for the clade containing Fare¬ 
well Spit and Waikuku, Houpoto, Opoutere, 
and Rarawa. 

There were no differences in the structure 
of the sclerites of the male pedipalpal bulb or 


the sclerites of the female external epigyne of 
the 22 museum specimens of L. katipo and L. 
atritus examined. 

DISCUSSION 

Although too few genetic samples were col¬ 
lected to gain a definitive view' on intra-spe¬ 
cific gene flow between L. katipo and L. atri¬ 
tus populations, an indication can be inferred 
from the results. The maximum interspecific 
pairwise distance among L. katipo and L. atri¬ 
tus populations was 0.97%, which was smaller 
than most pairwise distances found in the 
ND1 gene region between Nesticus spp. pop¬ 
ulations (Hedin 1997a). Although compari¬ 
sons between genera are not ideal, Nesticus is 
in a sister family to Latrodectus so some in¬ 
ference may be drawn. Low intraspecific pair¬ 
wise distances among L . katipo and L. atritus 
populations, therefore, indicate that popula¬ 
tions from which genetic samples were col¬ 
lected may not be genetically isolated or they 
have not yet undergone complete lineage sort¬ 
ing. 

Both maximum likelihood and parsimony 
analyses revealed that these taxa were para- 
phvletic. In addition to this, none of the 22 
museum specimens of L. katipo and L. atritus 
examined were found to differ in male and 
female genitalia. However; marked differences 
between L, katipo and L. atritus coloration 
and distribution (McCutcheon 1976; Forster & 
Kingsford 1983; Forster & Forster 1999) offer 
support for the current taxonomic designation 
of these species. Latrodectus are unusual 
amongst spiders in that their coloration ap¬ 
pears to be more useful than genitalia in sep¬ 
arating species (Levi 1983). It is possible that 
although L. katipo and L. atritus have not 
been observed mating, this does not preclude 
the possibility that these species may inter¬ 
breed. Moreover, if color variation between 
the species were related to an environmental 
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variable, such as temperature, differences in 
morphology and distribution may be ex¬ 
plained. However, we have only data from one 
mitochondrial gene region and the genetic dif¬ 
ferences between L. katipo and L. atritus may 
be due either to gene flow or incomplete lin¬ 
eage sorting. The validity of the species status 
of L. katipo and L. atritus could be further 
explored by the sequencing of more popula¬ 
tions with ND1, sequencing with other mito¬ 
chondrial genes (e.g., COI, which showed 
over three times more sequence divergence 
between L. katipo and L. hasselti than ND1: 
see Garb et al. 2004) or nuclear gene introns, 
the screening of microsatellites and/or detailed 
morphological examinations for other possible 
characters, especially including the mixed 
populations of the two species mentioned by 
Forster & Forster (1999). It would also be 
worth repeating the rearing experiments re¬ 
ported by Forster & Kingsford (1983) and 
Forster & Forster (1999) with large replicates. 
Until further work is undertaken L. katipo and 
L. atritus should be regarded as separate spe¬ 
cies. 

Forster (1995) postulated that the New Zea¬ 
land widow fauna had been genetically iso¬ 
lated from L. hasselti since the fragmentation 
of Gondwana 60-80 mya (Hayes & Ringis 
1973). If Forster’s (1995) hypothesis is true 
and L, hasselti and the New Zealand fauna 
have been isolated from one another for at 
least 60 my, the maximum pairwise distance 
(1.95%) between the L. hasselti and the New 
Zealand widow fauna suggests a rate of 
change the ND1 sequence of 0.0325% per 
million years, which is 70 times slower than 
mitochondrial sequence divergence reported 
in other arthropods (Brower 1994). Forster 
(1995) had also suggested that all Latrodectus 
species had a common theridiid ancestor be¬ 
fore the break up of Pangea 400 mya (Stevens 
1985), which would predate the earliest 
known Araneoidea fossil (Selden 1989) by 
270 my and the earliest spider fossil (Shear et 
al. 1989) by 20 my. Much of the present day 
distribution of Latrodectus is likely to be due 
to dispersal events (Garb et al. 2004) and the 
low genetic divergence between L. hasselti , L. 
katipo and L. atritus in this study and between 
L. hasselti and L. katipo in Garb et al. (2004) 
suggests that Latrodectus was not present on 
New Zealand when it separated from Gond¬ 
wana 60-80 mya. This assertion is supported 
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0.01 substitutions/site 

L katipo - Farewell Spit & Waikuku 


L. atritus - Rarawa 


L atritus - Opoutere 


L. atritus - Houpoto 


L. katipo - Kaitorete Spit 


L. katipo - Herbertville & Flat Point 


L. hasselti - Western Australia & Queensland 


L. he $ perns 


1.55 substitutions/site 

Figure 2.—Maximum likelihood tree, which has identical topology to one of the 12 equally parsimo¬ 
nious trees. Note the long branch of L. hesperus. There was not bootstrap support over 60% in the 
maximum likelihood analysis and none over 50% in the parsimony analysis. 



by growing evidence that the New Zealand 
spider fauna has been, and continues to be, 
influenced by the arrival of spiders from Aus¬ 
tralia (Vink & Sirvid 2000; Vink et al. 2002; 
Vink & Paterson 2003). Given that suitable L. 
katipo and L. atritus habitat has probably been 
present in New Zealand for a long time (Ste¬ 
vens et al. 1988) and that the genetic evidence 
indicates L. hasselti is a good disperser, it 
seems unlikely that L. katipo and L . atritus are 
recent arrivals to New Zealand. In the absence 
of datable fossil records, however, it is un¬ 
likely that the time L. katipo and L. atritus 
arrived in New Zealand will be precisely 
known. Overall, this lack of evidence for iso¬ 
lation between Australia and New Zealand 
since the Gondwanan break-up agrees with 
other recent studies of other New Zealand 
taxa, such as podocarp trees (Pole 1994), gal- 
axiid fish (Waters et al. 2000), hepialid moths 
(Brown et al. 1999) and various flightless in¬ 
sects (Trewick 2000). 

It is possible some gene flow may occur, or 


has occurred until recently between Australian 
and New Zealand Latrodectus populations. 
Greater intraspecific variation found among 
populations of L. katipo or L. atritus than be¬ 
tween the Australian specimens might have 
resulted from periods of glaciation or rising 
sea levels that restricted gene flow between L. 
katipo or L. atritus populations in New Zea¬ 
land, but are unlikely to have affected L. has¬ 
selti (Stevens 1985; Stevens et al. 1988; Nich¬ 
ols 2001; Trewick 2001). Moreover, Raven 
(1992) and Main (1992) suggested thatL. has¬ 
selti may have only recently been introduced 
to eastern Australia from South Australia, 
which would explain the lack of genetic var¬ 
iation between the L. hasselti specimens. 

Although the ND1 mitochondrial gene re¬ 
gion has previously been used to examine in¬ 
tra-specific variation between spider popula¬ 
tions (Hedin 1997a; Hedin 1997b; Masta 
2000; Johannesen et al. 2002; Maddison & 
Hedin 2003; Masta & Maddison 2002), this 
gene region did not evolve fast enough to pro- 
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vide the definition required to examine gene 
flow between populations of L. katipo , L. atri- 
tus or L. hasselti. Moreover, the low number 
of samples examined in this project also made 
it difficult to gain a definitive view of intra¬ 
specific gene flow. These problems might be 
overcome if sequence data from other faster 
evolving mitochondrial gene regions, such as 
COI (see Garb et al. 2004), or microsatellites 
were used and more samples were examined. 
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